Topological Defects Formation after Inflation on Lattice Simulation 
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We consider the formation of topological defects after inflation. In order to take into account 
the effects of the rescattering of fluctuations, we integrate the classical equation that describes 
the evolution of a complex scalar field on the two-dimensional lattice with a slab symmetry. The 
growth of fluctuations during preheating is found not to be enough for defect formation, and rather 
a long stage of the rescattering of fluctuations after preheating is necessary. We conclude that the 
topological defects are not formed if the breaking scale r\ is lager than ~ (2 - 3) X 10 16 GeV. 



PACS numbers: 98.80.Cq, 11.27. +d 



I. INTRODUCTION 

Inflation was invented in order to overcome sev- 
eral problems in the standard hot big bang universe, such 
as the flatness and horizon problems. Moreover, harmful 
topological defects created before inflation were diluted 
away, and are almost absent in the present universe. If 
some topological defects such as monopoles or domain 
walls are produced after inflation, they are disastrous 
since they will soon dominate the energy density of the 
universe. Therefore, the reheating temperature cannot 
be as high as the grand unified theory (GUT) scale so as 
not to produce the GUT monopoles. 

It was recently recognized that topological defects may 
be formed during preheating |3|-|| . Preheating is the very 
beginning of the reheating process, and occurs because 
of the parametric resonant effects . Its essence is an 
induced effect in the sense that the presence of the pro- 
duced particles stimulates further decay of the coherently 
oscillating scalar field into those particles jl(| . This phe- 
nomenon is thus peculiar to bosonic particles that obey 
the Bose-Einstein statistics. 

In the preheating stage, very large non-thermal fluc- 
tuations are produced, (S(j) 2 } ~ c 2 M 2 where M p is the 
Planck mass and c = 1CT 2 - 1CT 3 f|j]. These fluctu- 
ations change the shape of the effective potential of the 
field 4> to restore its symmetry if the potential V((j>) is of 
spontaneous symmetry-breaking type. Later, when the 
amplitude of these fluctuations is redshifted away by the 
cosmic expansion, the symmetry is spontaneously bro- 
ken and topological defects may be created. Thus, the 
mechanism for producing the topological defects is some- 
what similar to the Kibble mechanism in high tempera- 
ture theory. However, we showed in Ref. || that the 
amplitude of the fluctuations does not grow larger than 
that of the homogeneous mode for V(<j>) = A(</> 2 — ry 2 ) 2 , 
where A is the self-coupling constant. It implies that 
the symmetry restoration does not take place and hence 
the topological defects are not produced. Therefore, we 
considered that topological defects are formed through 
another mechanism which is based on the following two 
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facts. The first one is that there are two minima of the 
effective potential in the radial direction in the simple 
U(l) theory. The second is that some fluctuations ex- 
ist initially at the preheating stage. They are produced 
during inflation and stretched far beyond the horizon to 
become classical fields. Therefore, the initial condition 
for the homogeneous field in some region in the universe 
is different from that in another region in the universe. 
The dynamics of each region in the universe is thus dif- 
ferent and the final value of the field (the minimum into 
which the field settles down) is different, which leads to 
the defect formation. 

However, in Ref. j^j, we used the Hartree approxima- 
tion so that the interactions between particles with dif- 
ferent momenta are not fully included. As we will see 
later, rescatterings pd|dl5| are quite important to deter- 
mine the dynamics of the defect formation as suggested 
by Refs. (l^JT^]. To deal with the effect of the rescat- 
tering, we must go on to the lattice simulation [[ll] [l|. 
Therefore, in this paper, we integrated the equations de- 
scribing the dynamics in a real space on the lattice. There 
are some difficulties in the simulations on lattice. One of 
the difficulties is how we can implement the quantum 
fluctuations of the fields on the lattice. We rely on the 
idea of Ref. |ll[] and do not discuss in detail (actually, it is 
not sensitive to the initial conditions for the fluctuations 
if the order of magnitude is properly taken The 
other one is the box size and the lattice size. The lattice 
size must be small enough for identifying the defects. The 
box size should not be taken too much smaller than the 
horizon size at the end of the simulation, otherwise we 
miss the defects whose density is 0(10) per horizon. We 
also take the proper lattice and box sizes so that we can 
see the resonance effects. From these requirements, we 
calculate the evolution of a complex scalar field in a two- 
dimensional lattice abandoning to calculate it in three 
dimensions because of the lack of machine memories. 



II. MODEL 



Let us consider the wine bottle potential 
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FIG. 1. Evolution of the real part of the homogeneous 
field for V = 10 16 GeV. 



FIG. 2. Evolution of the fluctuations for rj = 10 16 GeV. 



(i) 



and 



where $ is the complex scalar field, A is the very small 
self coupling constant, and r\ is the breaking scale. This 
model has a global U(l) symmetry, and cosmic strings 
are formed when the symmetry is spontaneously broken. 
The equation of motion becomes 

$ + 3H$-47V 2 $ + A(|$| 2 -77 2 )<I> = 0, (2) 

or 

where H is the Hubble parameter and the dot denotes 
differentiation with respect to time t. Rescaling as 



i(r)dr = y/\$ a(0)dt, 
$o(r) 



*oo(0)' 
£ = \/A$ o a(0)a;, 

where $ = |$(0)|, then, setting a(0) = 1, Eq.(§) be- 



(3) 
(4) 
(5) 



comes 
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-^-V 2 ^+(M 2 -ryVV = 0, (6) 



where 77 = i]/§o and the prime denotes differentiation 
with respect to r. The second term of LHS can be omit- 
ted since the energy density of the universe behaves like 
radiation at the early time and also the scale factor be- 
comes very large later. Therefore, we assume that the 
universe is radiation dominated. In this case, the rescaled 
Hubble parameter, h(r) = H(t)/>/X^q, an d the scale 
factor a(r) become M 



Kr) = ^a- 2 (r), 



a(r) 



V2 



(8) 



(7) 



respectively, when $ is assumed to be an inflaton (even 
if $ is not an inflaton, results are the same as in the 
case of rescaling the breaking scale in an appropiate way, 
see Ref. ||) and we put a(0) = 1. Since the physical 
length and the horizon grow proportional to a and a 2 , 
respectively, the rescaled horizon grows proportional to a. 
The initial length of the horizon is 4(0) = 3/V2 w 2.12. 
Therefore, the box size should be larger than the horizon 
size at the end of calculation, lySj) — 4(0)a(r). 

The width of the topological defect is (\/A?7) _1 which 
corresponds to (^(r)) -1 in the rescaled scale. Since it 
decreases with time, one lattice length should be at least 
comparable with the defect width at the end of the cal- 
culation. 

Furthermore, in order to see the resonance effects, both 
the lattice size should be small enough and the box size 
should be large enough. Typical resonant momentum is 
k ~ n/A^o a t the beginning. Since, in the A|<&| 4 theory, 
the resonance band does not change because of cosmic 
expansion, the lattice size should be smaller than Ax ~ 
^res ~ (V^^o) -1 - It corresponds to ^ 1 in a rescaled 
unit. In addition, it is necessary to have the box size 
large enough for a good resolution in the resonance band. 

Three requirements lead us to take very large lattice 
sizes at the expense of one dimension in space. Thus we 
integrate Eq. (g) on a 4096 x 4096 lattice with a slab sym- 
metry in a ^-direction. Therefore, the defects that we see 
on the lattice are the cosmic strings stretched infinitely 
along the ^-direction. 

For the initial conditions we take 
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FIG. 3. Spectrum of the fluctuations for 77 = 10 16 GeV. 
These lines shows r = 800, 500, 200, 130, 100 from top to bot- 
tom. 



ip x = Reip(0) = 1 + Sip x (x), 
ip y = Imyj(O) = Sip y (x), 



(9) 
(10) 



where the homogeneous part comes from its definition 
(we call x-direction for real direction and y for imag- 
inary), and 8ip XtV {x) is a small random variable rep- 
resenting the fluctuations (we owe the fact that the 
semiclassical description can be used to Ref. jllj). We 
also put small random values for velocities. We inte- 
grated Eq.(^|) until the time when the horizon grew as 
large as the box size for five different breaking scales, 
j] = 10 16 ,2xl0 16 ,3xl0 16 ,5xl0 16 , and 10 17 GeV. For ex- 
ample, we integrate until r max — 2000 for rj = 10 16 GeV. 
To check the accuracy of our code, we calculate the en- 
ergy density of the field for the V(&) — A|$| 4 model, and 
it conserves within 0.36% during preheating and 0.07% 
after. (Notice that the energy density is not constant in 
the wine bottle type potential.) We identify the defect if 
there is a winding number at that point. 



III. RESULTS 

We will first show the general feature of the dynam- 
ics taking 77 = 10 16 GeV as an example. We take the 
lattice size A£ = 0.5 and (^(x)) 1 / 2 = 10~ 7 - 10~ 10 
in this case. Figure [l] shows the evolution of the real 
part of the homogeneous field. Preheating takes place 
at t ~ 100 and ends at t ~ 150. As can be seen in 
Fig. ^], the fluctuation in y-direction grows first and the 
x-direction follows it. This is because the field feels no 
potential in the phase direction, which means that the 
y-direction corresponds to the Nambu-Goldstone mode. 



FIG. 4. Evolution of the fluctuations for rj = 3 x 10 16 GeV. 

During preheating, only those modes in the resonance 
band grow exponentially through the parametric reso- 
nance effect (see Fig. ||). Since the mode equations of 
the fluctuations in x- and y-directions can be written as 



64 + [kl 
Sy'k + [k 



3x 2 



rj 2 a 2 }5x k 
■ rj 2 a 2 ]5y k 



0, 
0, 



(11) 
(12) 



where back reactions are neglected and k = fc/VA^O: the 
resonance mode for each direction is 3/2 < k 2 < \/3 and 
k 2 < 1/2. As we mentioned above, the fluctuation in 
y-direction grows faster. In the actual calculation, the 
fastest growing mode is found to be k ~ 0.47, which co- 
incides with the initial horizon size, ~ £h(0). We also 
change the initial amplitude of fluctuations. The differ- 
ence only appears in the time at the end of preheating 
T p h but its change does not shorten the duration of the 
rescattering so much [see Eq.([l3|) below]. 

After the exponential growth of the fluctuations (the 
preheating stage), there is a rather long stage of oscil- 
lations of the homogeneous mode with small amplitude, 
which ends when the field settles down into the minimum 
of the potential at r ~ 700. During this stage, rescatter- 
ings of the fluctuations become important. We can see 
that the amplitude of the fluctuations slightly grows in 
Fig. ^|, but redistribution of the spectrum is much more 
efficient. 

Let us consider the criterion for the defect formation. 
The most simple idea is to compare the time when the 
amplitude of fluctuations grows because of the paramet- 
ric resonance (r p h) with that when the field settles down 
into the minimum of its potential (r/ /;). More pre- 
cisely, Tph is the time when the homogeneous field de- 
cays through the parametric resonance until the back re- 
actions of the created fluctuations cannot be neglected. 
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FIG. 5. Spectrum of fluctuations for r\ = 3 x 10 16 GeV. 
These lines shows r = 250, 200, 150, 100 from top to bottom. 



This usually happens when the amplitude of fluctuations 
becomes as large as that of the homogeneous mode. It is 
estimated from 



Stp(0) exp(nT ph ) ~ 1, 



(13) 



where fj, is the effective growth exponent. In the case 
of Eq.(p"2|), the maximum value of /! is \i m ax — 0.147 
CSffa. Since T p h — 135 for 5(p(0) = 10~ 7 which we take 
in Fig. |, T ph ~ 150 for 5<p(0) = 10" 8 , and T ph ~ 180 for 
5<p(0) = 10~ 10 in our calculations, we get [i e ff ~ 0.12. [j] 
On the other hand, the amplitude of the homogeneous 
field decreases due to cosmic expansion, and at r = Tf a u , 



o(0) 



$0 



(14) 



Thus Tph ^ Tf a u for 7] > 6.7 x 10 16 GeV. In this case, 
preheating ends before it is fully developed. Therefore, 
topological defects are not formed, in this criterion, if 
77 > 6.7 x 10 16 GeV. In fact, we find no topological defects 
for 77 = 10 17 GeV. However, this is not a good criterion 
as we will see below. 

Actually, the amplitude of the fluctuations becomes 
a little smaller than that of the homogeneous mode at 



^Heff is estimated from Eq,([l3|). Since the resonance does 
not occur from the very begining and becomes effective a 
little later, the more conventional estimation is fj, ~ 0.141. 
Actually, the square of amplitude of the fluctuation grows 
from ~ 3 x 10 -15 to ~ 0.1 during resonance (At = 110) for 
5<fi(0) = O(10~ 7 ). 



FIG. 6. Evolution of defects for 77 = 2 x 10 1B GeV on lattice 
at (a) t — 1000 and (b) r = 1200. Dots and circles denote 
defects and anti-defects, respectively. 



most even if the breaking scale is lower (x 2 ~ 0.25 while 
(S(p 2 ) ~ 0.2), and they cannot affect the dynamics of 
the field wholly. This fact was also seen in the previ- 
ous work using the Hartree approximations which con- 
tain the forward scatterings only J^]. Therefore, we con- 
clude that the symmetry cannot be restored only through 
the parametric resonance. However, after the preheating 
ends, the mode-mixing (rescattering) of the fluctuations 
becomes important for determining whether or not the 
topological defects are produced. 

Now we are going further into the detail of the rescat- 
tering. At the end of preheating, the occupation numbers 
of the fluctuations are very large only in the resonance 
band. The typical time scale for the rescattering is es- 
timated as follows |14|. The cross section for a scatter- 
ing is naively given by o\ ~ \ 2 /E 2 . Here E ~ is 
the typical energy (momentum) of the produced particle. 
However, since a huge amount of particles are present af- 
ter preheating, the cross section is much larger because 
of the Bose enhancement. Therefore the phase space 
distribution function (the number density in /c-mode) 



,/E 3 should be multiplied to o\. After all, the 



scattering rate becomes 



(15) 



Using Eqs.(||) and (g), we get a typical time scale for the 
rescattering T sca t ~ 0(1). Then the occupation number 
at higher momenta will be significant within At ~ 30. 
Notice that the effective rescattering is essentially the 
same mechanism as the parametric resonance in the sense 
that these are stimulated effects. The process 5(p(k res ) + 
6ip(k res ) — > Stp(ki ow ) + Stp(khigh) nrs t occurs. Typically, 
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FIG. 7. Phase of the field around the string and anti-string 
at x ~ 1510 and y ~ 3830 in Fig. |(d). 

kiow — because of the Bose enhancement, since the 
occupation number with k = (homogeneous mode) is 
large. In this picture some major processes are as follows: 

8ip(k res ) + 5ip(k res ) — > ip + 5(p(2k res ), 

8<f(2k res ) + ip — > 8(p(k res ) + <fyj(fe res ), 

8<p(2k res ) + dip(2k res ) — > 8tp(k res ) + Sip(3k res ), 
8<p(3k res ) + ip — > 8tp(k re s) + <fyj(2fc res )- (16) 

These are efficient because the produced fluctuations of 
mode k (the first terms of the RHS's) have already large 
occupation numbers. However, these processes are the 
zero-sum game. If high momentum fluctuations are pro- 
duced, the homogeneous mode is also created and cannot 
decay further. 

There are also processes which do not feel the bose 
enhancement in the first place such as: 

5ip(k res ) + ipo — > 8ip(ki) + 8ip(k res - fci), 

5<p(2k res ) + ip —> 8tp(k 2 ) + 8ip(2k res - k 2 ). (17) 

The field is knocked off from the homogeneous mode 
because of these processes which make the spectrum 
smooth. Combining the above two types of rescatter- 
ings, the occupation numbers of the mode outside the 
resonance band grow large enough for the zero mode to 
decay into all the modes. The decay of the homogeneous 
mode completes in such a complicated way, and the topo- 
logical defects will be formed through the Kibble-like 
mechanism later. Therefore, the rescattering timescale 
for the defect formation is considerably longer than that 
given in Eq.(|l5|). 

Now let us go on to the results of the rescattering ef- 
fects. As we can see in Fig. ||[ the fluctuations in the 



FIG. 8. Changes of the number of the defects for 
r\ = 10 16 GeV with various initial conditions of fluctuations. 
The solid line denotes the average. 

resonance band develops in the first place, and later it 
moves to both higher and lower momentum modes by 
the rescattering. The shape of the spectrum becomes 
very broad and smooth, and further decays of the homo- 
geneous mode into all the modes can be seen because 
all the amplitudes grow large. We show the case for 
r] = 3 x 10 16 GeV in Figs.g and || Although the 
amplitude of fluctuations becomes as large as that of 
T] = 10 16 GeV, both higher modes and lower modes do 
not develop so much in comparison with Fig. ||. There- 
fore, it is not enough time for the rescattering to produce 
the topological defects when rj — 3 x 10 16 GeV. Actually, 
we find no topological defects at the end of the calcula- 
tion for this breaking scale. Consequently, it is necessary 
to have a somewhat long stage of the rescattering in or- 
der to produce the defects. We find that the topological 
defects forms for 77 < 2 x 10 16 GeV at most and can thus 
conclude that duration of rescattering needs t > 200. 

Here we investigate cases for 77 = 10 16 GeV and 77 = 
2 x 10 16 GeV in some extent. In the latter case, there are 
only a few strings in the horizon size (the whole region 
at the end of calculation), as seen in Fig. |(]. The identifi- 
cation of the defects is done by observing the winding at 
the space point on the lattice as is seen in Fig. [7| Since 
half of strings are in fact anti-strings and each pair of the 
string and anti-string is very close to each other, it will 
annihilate and disappear very soon, which means that the 
case of 77 = 2 x 10 16 GeV may be harmless for cosmolog- 
ical history. Actually, the numbers of defects at the end 
of calculations are 2, 4, 2, and 16 for four runs. On the 
other hand, for the case of 77 = 10 16 GeV, even though 
the number of strings in the horizon decreases as time 
goes on, they still exist in a certain amount (see Fig. ||). 
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FIG. 9. Evolution of the defects for r) = 10 16 GeV on the 
lattice at (a) r = 1300, (b) r = 1600, (c) r = 1800, and (d) 
r = 2000. Dots and circles denote defects and anti-defects, 
respectively. 

We expect that it will settle down to some kind of a scal- 
ing solution later. The evidence that topological defects 
will survive annihilations is that the average number of 
strings at the end time of calculation (r = 2000) is a 
little larger and the positions of strings and anti-strings 
are not always pair-like, which is seen in Fig. |^ (d) as 
opposed for 77 = 2 x 10 16 GeV in Fig. ||. 

IV. CONCLUSIONS 

We have considered the formation of topological de- 
fects after inflation. To this end, we have integrated 
the equation which describes the evolution of a complex 
scalar field with the potential V(4>) = A(|$| 2 - rj 2 ) 2 /2 
on the two-dimensional lattice with a slab symmetry in 
the z-direction in order to include both parametric reso- 
nance effects and rescatterings of fluctuations. We have 
found that fluctuations produced during preheating do 
not lead to the defect formation even if their amplitudes 
become the same order of magnitude as that of the ho- 
mogeneous mode at the end of preheating. We also have 
found that the effects of the rescattering of the fluctu- 
ations are essential and a rather long duration of the 
rescattering is necessary for the topological defects to be 
formed as suggested by Refs. [Uj,[17|. Fluctuations in the 
resonance mode scatter off each other to make fluctua- 
tions with higher and lower momenta. They also knock 
the field off from the zero mode to broaden and smooth 
the spectrum, which in turn stimulates further decays of 
the homogeneous mode into all the modes. Therefore, 
the topological defects are formed in a similar way as the 



FIG. 10. Evolution of the fluctuations for 77 = 10 16 GeV 
using the Hartree approximation in two dimensions in space. 

Kibble mechanism. These processes take a lot of time so 
that it is necessary for the field not to fall into the min- 
imum of the potential too soon after preheating. This 
constrains the value of the breaking scale. From our cal- 
culations, we conclude that the topological defects are 
not formed if 77 > (2 - 3) x 10 16 GeV. 

The constraint on the breaking scale is higher than 
that of our previous result in Ref. || where the Hartree 
approximation is adopted. In both cases, the amplitude 
of fluctuations at the end of preheating is not larger 
than that of the homogeneous mode (compare Figs. ^ 
with . There are some differences between the calcu- 
lations with the Hartree approximation and the lattice 
simulation during preheating. The fluctuations in the 
y-direction have a similar shape, but the amplitude ob- 
tained with the use of the Hatree approximation is a little 
smaller. The fluctuation in x-direction in the calculation 
with the Hatree approximation does not grow until a very 
late time, which differs completely from the lattice one. 
The reason is that most of the rescattering effects are 
neglected in the Hartree approximation while the rescat- 
tering is very efficient for producing fluctuations in both 
x- and y-directions in the lattice calculations. Therefore, 
it is not surprising that the critical breaking scale for the 
formation of topological defects becomes higher in the 
lattice simulations. 

Note added. After submission of this paper, we no- 
ticed a paper, Ref. JTs| ] , in which numerical calculations 
were done in three dimensions, and similar results were 
reported. We agree with the authors of Ref. [^9| that 
rescatterings are less effective in two dimensions to some 
extent. However, there is little difference in the resonance 
effects. Actually, /1 ~ 0.141 in ours [see the footnote be- 
low Eq.@] while ~ 0.147 in Ref. jT||. Moreover, the 
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maximum value of fluctuations is almost the same. The 
advantage of two-dimensional simulations is that the box 
size can be taken as large as the horizon size at the end 
of calculations while the lattice size can be taken small 
enough to identify defects. Therefore, we can estimate 
the average numbers of defects per horizon size. 
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